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ABSTRACT 

We develop a new analytic solution for the steady-state structure of a thin accretion 
disc under the influence of a magnetic field that is anchored to the central star. The 
solution takes a form similar to that of Shakura and Sunyaev and tends to their solution 
as the magnetic moment of the star tends to zero. As well as the Kramer's law case, 
we obtain a solution for a general opacity. The effects of varying the mass transfer 
rate, spin period and magnetic field of the star as well as the opacity model applied 
to the disc are explored for a range of objects. The solution depends on the position 
of the magnetic truncation radius. We propose a new approach for the identification 
of the truncation radius and present an analytic expression for its position. 

Key words: accretion, accretion discs - stars: magnetic fields - stars: pre-main- 
sequence - stars: dwarf novae - stars: neutron 



1 INTRODUCTION 

There is a well known one dimensional solution for the 
structure of a thin accretion disc in the steady state. 
This solution, known as the Shakura-Sunyaev disc solu- 
tion (Sha kura fc Sunvaevll973h . consists of a series of seven 
equations. The full solution contains expressions for disc 
surface density, scaleheight, density, temperature, opacity 
and viscosity. These expressions depend upon the Shakura- 
Sunyaev alpha parameter, a dimensionless parametrization 
' of viscosity. In addition to the alpha parameter the solu- 
tion is expressed as a function of radius, mass transfer rate, 
stellar mass and stellar radius. 

We wish to reformulate the Shakura-Sunyaev disc so- 
lution to incorporate the effect of a torque from a stel- 
lar magnetic field. A solution for the structure of an ac- 
cretion disc in a magnetic field is desirable because accre- 
tion discs are frequently found about magnetic stars. Such 
magnetic fields can allow the transfer of angular momen- 
tum between the disc and star and can therefore signifi- 
cantly affect the structure of the disc as well as the spin 
evolution of the star (e.g . iBrandenburg fc Campbell Il998t 

lArmitaee fc Clarke! Il99 fjl. These effects can be important 
in yo ung stars (e.g. lKonigJll99lt |Anntege^Clar ke fc Toutl 

I1999T) . cataclysmic variables (e.g. lsTrlerikeretHuTl20o!3) and 
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X-ray binaries (e.g. iPringle fc Reesl|l972l : iRomanova et all 

2003). Accretion discs can exhibit outbursts due to the 
thermal viscous disc instability, and the interoutburst time, 
known as the recurrence time, can be increased by sev- 
eral orders of m agnitude by truncating the inner disc (e.g. 
iMatthews et aljl2004) . Accretion disc truncation by various 
mechanisms has also been invok ed to explain the ultraviolet 
(UV) lag in dwa rf novae (e.g. iMever fc Mever-Hofmeisteil 
Il994l : lKindll997l l. In this paper we treat magnetically en- 
hanced accretion, but not the magnetic propeller effect. 

A solu tion for the structure of a d i sc in a magnetic field 
is given bv IBrandenburg fc CamobelJ 1^998). By neglecting 
some mass transfer terms, they are able to find a detailed 
solution for the structure of the inner disc where a specific 
opacity model is assumed. In their solution magnetic pres- 
sure is considered in addition to magnetic torque, which im- 
proves the accuracy of the model for the inner disc regions. 
However, the azimuthal magnetic field is parametrized with 
respect to the shear between the disc and the stellar field us- 
ing essentially the same form as is used in this paper. Since 
the effect of mass transfer into the disc is partially neglected, 
the solution is most valid in the very inner parts of the disc, 
where the magnetic field dominates the disc structure com- 
pletely. A complementary solution is presented in this paper 
which, though less accurate in the centre of the disc, is con- 
tinuous for the entire disc and can be applied to a range 
of opacity models. Furthermore the solution derived below 
tends to the Shakura-Sunyaev solution in the case where the 
stellar magnetic field vanishes, and also at large radii. 
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The paper begins with a review of the necessary as- 
sumptions for the new solution. Then an approximate for- 
mulation for the magnetic torque is derived. This is applied, 
in conjunction with a new estimate of the magnetic trunca- 
tion radius, to the derivation of a disc model. This model is 
generalised in order that it is compatible with any opacity 
which takes the form of a power law in density and tem- 
perature. We then explore how variables such as the mass 
transfer rate, magnetic field and stellar spin influence the 
disc structure. 



2 THE SOLUTION 
2.1 Assumptions 

The derivation requires eight equations which follow in 
a form which is appropriate for our new magnetic case. 
The equations used here ar e similar to those presented by 
iFrank. King fc Rainel l)2002h in their derivation of the stan- 
dard Shakura-Sunyaev solution. First, density in the disc 
can be approximated by 

'=§' W 
where E represents surface density and H is the scaleheight 
of the disc. For a thin disc the scaleheight can be approxi- 
mated as 



H 



S R 3/2 



(2) 



(GM*) 1/2 

where c s is the local sound speed, R is the radial distance 
from the centre of the star in the plane of the disc and 
is the mass of the primary star. The universal gravitational 
constant is represented by G. This approximation is likely 
to break down close to the star, where a strong magnetic 
field may cause gas to stream along magnetic field lines. An 
ideal gas has an equation of state of the form 



2 P 
C s = — 



(3) 



in which P represents pressure. Pressure is composed of 
three parts in the magnetic case: gas pressure, radiation 
pressure and magnetic pressure. We shall neglect the con- 
tributions of radiation and magnetic pressure, as these are 
much smaller than gas pressure for most cases. This yields 

pfeT c 



P = 



fj,m u 



(i) 



where k is the Boltzmann constant, T c is the temperature in 
the mid-plane of the disc, the atomic mass unit is given by 
m u , and n is the mean molecular weight in the disc. This is 
frequently taken to be \i ~ 0.6. Conservation of energy gives 
jFrank et al.ll2002l) 

4aT c 4 9 ^GM„ ... 
~TT = 8" S -R^ ' (5) 
if it is assumed that the disc is Keplerian until the boundary 
layer and that the disc has a large optical depth. Here a rep- 
resents the Stefan-Boltzmann constant and r is the opacity 
of the disc material. At this stage the form of the function 
zvE is unknown, so that the equation must be kept in this 
general form. It is initially assumed that opacity may be 
approximated by Kramer's law such that 



T = K EpT c 



-7/2 



with 



5 x 10 24 cm 5 g- 2 K 7/2 



(6) 



Kramer's opacity is not a suitable model for all discs. Alter- 
native opacity laws are discussed in Section 12.51 Following 
the same approach as for the normal Shakura-Sunyaev so- 
lution, we assume that viscosity can be parametrized in the 
form 



v = aCsH 



(7) 



In addition, an expression for uT, must be obtained. In the 
classical, non-magnetic Shakura-Sunyaev solution the rela- 
tion 

1/2] 

(8) 



v M 
3tt 



1 



(!) 



holds, where J?* denotes the radius of the star. However in 
the magnetic case the relation is more complex. 

2.2 Magnetic torque 

In order to obtain an expression for i/E in the magnetic case, 
we integrate the continuity equation in the direction normal 
to the plane of the disc, which gives 

dt R OR W 
Similarly the radial component of the Navier-Stokes equa- 
tion becomes 



/ dv R dv R v 2 \ dp 
E {-dT +V *^R-R ] = 



,GM 4 



<:) 



3i? 3 / 2 OR 



R 3/2 uT. 



OR 



dR R? 

2 d(R 2 vT,v n ) 



(10) 



3R 3 

and the azimuthal component becomes 

dv. 



dR 



,_ Mr d(Rv v ) 

dt R dR 



R 2 dR 



R 3 uE 8 



dRKRJ 



(11) 



/a/ di\ 

^{dl +VR dR) 



+ EA 



(12) 



where v v and i>r represent the azimuthal and radial compo- 
nents of the velocity, respectively, and v is the viscosity in 
the disc. 

If we multiply equation Ijllfl by R, then we obtain an 
equation for the specific angular momentum, 

RdR [ R ^ 

where we have added an external torque acting on the 
disc, with A representing the injection rate of specific 
angular momentum and I indicating the current specific 
angul ar momentum. Th is follows the same procedure as 
iLin fc Pap algizffj 1 1986) , while a similar method is also em- 
ployed by Pringle (1991). 

Assuming that the disc is sufficiently cold and therefore 
that the dynamical time-scale is much smaller than the vis- 
cous time-scale, we can adopt a Keplerian approximation for 
the azimuthal motion such that v v = {GM* / R) 1 ^ 2 . Equa- 
tion il 1 2t can then be solved for the radial velocity in the 
disc, 



I'R 



d(R^vZ) +2A _R 1/2 



RV 2 E dR 



(13) 



^/GMl 

Inserting equation il l Ml into the continuity equation @ 
yields an evolution equation for the surface density of the 
disc, 
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<9E _ 3__d_ 

dt RdR 



R 



1/2 



d(R 1/2 vZ) 



dR 



1JL 
RdR 



2AE 



R^ 2 
VGAL 



(14) 



where the right hand side is composed of a diffusion term, 
and an advection term that is due to the external torque. 
To parametrize the specific torque A, we can write 

A = i = ^, (15) 

t\ tA 

where tA is the time-scale on which the local disc material 
gains angular momentum. 

In the present case, the source of the torque is the mag- 
netic interaction of a rotating, magnetic star with a partially 
ionised disc. The torque time-scale is therefore equivalent to 
the magnetic time-scale, tA ~ imag- 

The inclusion of a magnetic field in such a hydrody- 
namic system introduces two additional terms to the Euler 
equation: a magnetic p ressure term and a magnetic tension 
term, e.g. lDendvl (Il990h . The magnetic pressure term is neg- 
ligible where B is small. The magnetic tension term can be 
expressed as 



pr c V 



B z B<f 
4ir 



(16) 



where r c represents the local radius of curvature of the 
field lines as a result of the torque, while B z and B v rep- 
resent the vertical and azimuthal components of the mag- 
netic field respectively. We us e the approximation r c ~ H 
^Pearson. Wvnn fc Kingll997t) . The ratio of vertical and az- 
imuthal field strengths is related to the shear between the 
disc and the magnetic field. If it is assumed that the field 
rotates wi th the star then this r atio can be expressed in the 
form (e.g. iLivio fc Pringielll992h 



Bp 



(»k - op 



(17) 



where Q k represents the angular frequency of a body in a 
Keplerian orbit at a given radius and f2* = Ob denotes the 
angular frequency of the star and therefore of the magnetic 
field at all radii. The magnetic time-scale can now be defined 
in terms of magnetic acceleration and the Keplerian velocity 
by the relation 



tn 



Rn k 



4ttR 2 P H 



B'i 



(fik - 00 



(18) 



Finally, the volume density can be related to the surface 
density using equation Q and for a dipole field, we have 
B z ~ |/x| R~ 3 where fi is the magnetic moment of the star. 
These relations lead to the following expression for the mag- 
netic time-scale in the disc: 



tn 



AtyEVGM:R 11/2 



([£ 



3/2 



(19) 



where R co represents the corotation radius, at which a body 
in a Keplerian orbit will revolve about the star at the same 
angular frequency as that with which the star spins. If the 
field is moving more rapidly than the disc, which occurs 
outside the corotation radius, then the disc gains angular 
momentum and is pushed out to greater radii. Inside the 
corotation radius the disc moves more rapidly than the field 
and hence loses angular momentum and is accreted more 
rapidly. These regions are known as magnetic propeller and 



magnetic accretion regimes respectively. Because the mag- 
netic propeller, at least in this one-dimensional treatment, 
completely suppresses accretion, the disc cannot reach a 
steady state in the propeller case. Our steady-state solu- 
tion must therefore be confined to the magnetic accretion 
case. The re sult in equat ion l|19H is almost identical to that 
adopted by ILivio fc Pringld II1992T) . It is possible to use 
other prescriptions for the m agnetic time-scale, for exam - 
ple the dia magnct i c case /e.g.lWvnn. King fc Horn "el ll997h . 
ICampbell fc Heptins tall ( 1998a b) discuss the effects of dif- 
ferent assumptions for diffusivity. In general however t mag 
can usually be parametrized in the form 



2E 



^(7+2) 







([it:. 



3/2 



(20) 



where 7 and /3 are parameters determined by the magnetic 
interaction model. In the fully magnetised case represented 
by equation 1191 we have 7 = 7/2 and is determined by 
the mass and magnetic moment of the star according to 



(21) 



In fact the precise value of /3 may vary by a factor of order 
unity. This is because the choice of a fully magnetised disc 
as well as the adoption of the relation r c ~ H are both 
simplifying assumptions. 

By considering equation (1201 . a general expression for 
the advection term in equation 1141 is obtain ed which is 
similar to that found bv lLivio fc Pri ngle ( 1992) 



RdR 



2AE 



R?l 2 



VgmI 



p d 




r r ~ 


3/2 \ 1 


RdR 


R-i \ 


-Rco- 



.(22) 



Equation 1141 can now be solved analytically for the steady 
state. By definition, in the steady state, the mass transfer 
rate M — —2-kR'Svb, is a constant throughout the disc. By 
setting the left hand side of equation 1141 to zero and in- 
tegrating the remaining terms, this condition can be used 
to derive a general solution for i^E in the steady state such 
that, 



3^ + 3(2-7) 



R 



Rc- 



3/2 



(2-7) 
(2- -7) 



+CR- 



,(23) 



where C is an arbitrary constant. 



2.3 Truncation radius 



The truncation radius R t is the radius within which the ac- 
cretion rate is greatly enhanced by the magnetic field, and 
the surface density falls rapidly to a value close to zero. 
The position of this radius is of particular interest as it 
can provide a boundary condition for equation 1231 . Math- 
ematically the truncation radius can be identified with the 
largest real root of equation (I2.'!l . In the standard Shakura- 
Sunyaev solution the boundary condition vYj — is applied 
at the boundary layer, near R = 7?*, since there is no vis- 
cous transport at that point. Physically, in the magnetically 
truncated case, the surface density at radii smaller than the 
truncation radius must always be zero. Since equation il-Ml 
does not take this form, it is clear that the equation becomes 
unphysical where R < Rt- It is therefore difficult to justify 
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Tension truncation radius 
Alfven (pressure truncation) radi 




Figure 1. Plot to show how the Alfven and tension truncation radii, as measured in units of corotation radii, vary as functions of the 
parameter Q, denned in equation 1261 . The tension truncation radius is plotted for 7 = 7/2. Note that the tension truncation radius 
never exceeds corotation. 



the use of the usual boundary condition in this case, since it 
will lie in an unphysical regime. If a truncation radius exists 
and can be identified, then it follows that vYj — at R — R t 
may be used as a boundary condition. 

There have been several attempts to locate the trunca- 
tion radius due to the stellar magnetic field. The simplest 
of these estimates is the Alfven radius. This is the radius at 
w hich magnetic press ure is equal to ram pressure. According 
to lFrank et alJ (120021) the Alfven radius can be expressed by 



Rah 



5.1 x 10 s M" 2/7 Mr 



1/7 4/7 
' ^30 ' cm 



(24) 



where Mie represents the mass transfer rate in units of 
10 16 g s _1 , Mi is the stellar mass in solar masses and /X30 
is the stellar magnetic moment in units of 10 30 G cm 3 . This 
equation is not sensitive to the spin of the star. Indeed by 
inspection we see that the Alfven radius can be related to 
P/M so that 

Rah = 1-25 (/3/A/) 2/7 cm . (25) 

A more sophisticated approach than that of the Alfve n 
radius is adopted by iBrandenburg fc CampbelH lll99Sl) . 
who include magnetic tension and magnetic pressure in 
their calculations. Since some mass transfer terms are ne- 
glected, their solution is most v alid very close to the star. 
IBrandenburg fc Campbell! (jl99fil) do not provide an analytic 
formulation for the truncation radius, but according to their 
data table the radii that they calculate are all close to the 
Alfven radius. 

The approach adopted here, which is consistent with 
section |2~^1 in that it treats the magnetic tension only, is to 
solve equation l)2Mfl for R t by taking advantage of the fact 
that the inner accretion disc is empty in the case of mag- 
netically enhanced accretion. If E (Rt) = 0, and if annuli at 
smaller radii also have a zero surface density then, in addi- 
tion, i/S = at R = Rt and at immediately smaller radii. 
It is reasonable to conclude from this that d (vll) / dR = 
at R = Rt. Physically this can be interpreted as the ef- 
fect of a very rapid magnetic accretion. The inner edge of 
the disc loses angular momentum on a short time-scale and 



is dragged towards the stellar surface so rapidly that the 
inner disc has a negligible surface density. These two condi- 
tions for uY, are sufficient to obtain a relationship between 
mass transfer rate, /3, truncation radius and corotation ra- 
dius from equation 1231 so that 



M 
J 



2-kR: 



(—) 



3/2 



(26) 



This result can already tell us something about the be- 
haviour of the truncation radius in this approach. It is clear 
that when M//3 = 0, and the magnetic effect is at its most 
extreme, we have Tit = Rco- It is not possible for R t to 
exceed R co providing that M//3 > 0. This agrees with our 
previous assertion that the magnetic propeller does not have 
a steady-state solution. Of course when P/M = there is 
no truncation and we find that R t — 0. In reality however, 
since the inner radius of the disc cannot be smaller than 
the stellar radius i?*, the minimum value for the truncation 
radius is always R t = R*. It is also clear that, according to 
this formulation, Rt is independent of the disc opacity, and 
of a except in as much as M is related to a. There is no 
general analytic solution for Rt from equation , however 
it is possible to obtain a solution for certain values of 7. In 
order to do this, and to interpret the result, it is useful to 
make a parametrization of equation (I2(it . We define 

l_ 

-Mj 



= 2tt.R: 



(27) 



Combining this with equation l|26|l yields 



1 



' Rt 



,3/2 



(28) 



The tension truncation radius given by equation 1281 . and 
with 7 = 7/2, is plotted in Fig.Q The Alfven radius, which 
is also plotted, is smaller than the tension truncation ra- 
dius for much of the regime of interest. Where the Alfven 
radius is greater than the tension truncation radius the mag- 
netic pressure may be said to dominate. However, when the 
Alfen radius exceeds the corotation radius, it cannot repre- 
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sent the truncation radius of a steady-state accretion disc. 
For this reason the truncation radius given by equation 128H 
is preferable. 

The great advantage of the Q parametrization is that 
once the value of Q has been established from equation 127H 
it is straightforward to read off the value of Rt/R co , which 
can only vary from zero to unity. This can be done either 
from Fig. 0or from similar plots for different values of 7. 
This means that only one plot is required for all the param- 
eter space of interest. For the fully magnetised case, where 
7 = 7/2, an exact analytic solution cannot be found for 
Rt/Rco from equation In the region where Rt/R co is 

close to unity, the 7 = 7/2 case can be closely approxi- 
mated by setting 7 = 3, which yields a quadratic equation 
in (Rt/Rco) 3 ^ 2 ■ This, of course, has an exact analytic solu- 
tion, and since this is very close to the solution for 7 = 7/2 
for all of the region of interest, we can make for that case 
the approximation 



Rt 

Rco 



f 1 4/^+1-1 



2/3 



(29) 



While the fractional error in this approximation becomes 
larger when Rt/R co approaches zero, the absolute error is 
always less than 0.05, and, in any case Rt/R co can never 
become very small as physically it should be limited by R*. 

2.4 Synthesis 

The re sult in equation il'i.'it was obtained in lMatthews et all 



(2004). However in that work the standard non- magnetic 
Shakura-Sunyaev viscosity prescription was substituted into 
the equivalent of equation (IL'MI in order to obtain a surface 
density profile. This is clearly not self-consistent. In this pa- 
per no assumption of the form of v wil l be m ade other than 
that in equation (J. iMatthews et alJ (|2004h also used an 
arbitrary inner boundary condition such that dT,/dR = 
at R — 7?*. In the following derivation the inner boundary 
condition E (Rt) = 0, which is equivalent to the condition 
z/E = at R — R t , will be used, as discussed in Section f2. 31 
The value of Rt can be obtained using equation 1291 or by 
another method. In any case, where R t < R+, the boundary 
condition E (i?* ) =0 should be used instead, i n agreement 
with the usual approach fe.g. lFrank et al]l2002t) . The above 
boundary condition for R t > i?* yields an expression for the 
constant C, when applied to equation 1231 . so that 



i/E : 
where 
h = 



M 
3tt 



1 - 





l/2\ 


.R. 





f3R~ 



3 (7 - 2) 



h = 



M 4 
in 



(30) 



tion © as expected. Substituting relation 1301 into equation 
gives 



4aT c ' 



9 GMi, M 



,.4 

/ v 



(34) 



3r 8 R? 3tt" 

Following a procedure similar to that of lFrank et alJ l|2u 
we can now obtain a full analytic solution for the thin ac- 
cretion disc, to which a magnetic torque is applied. Solving 
the set of algebraic equations lHJ-Q, ©,10, lj30|l and (Sj} 
yields the following expression for surface density: 



/j32o_\ 1/10 /^m„\ 
V27k / \ k ) 



3/4 



-4/5 ( 



\ R 3 



1/4 



\ 3tt 



7/10 



(35) 



If n = 0.6 then the result can be rearranged into a more 
familiar form so that 



3/4,14/5 7/10 -2 
' f ' v ' g cm 



(36) 



E = Z.7a- A/6 MU w Ml ,A R-*i*f 

The mass transfer rate in units of 10 16 g s _1 is represented 
by Mi6 and Mi is the mass of the accreting star in solar 
masses. The radial distance from the star, in the plane of the 
disc, is represented by R in units of 10 10 cm. The numerical 
coefficient is sensitive to the values adopted for both /x and 
for opacity. Other than this the result in equation 136H differs 
from the usual form only by the presence of the correction 
term v, and the substitution of R t for throughout. It is 
instructive that the correction term differs from unity by an 
amount which is proportional to the ratio of /? to the mass 
transfer rate. This means that, as would be expected, a suf- 
ficiently high mass transfer rate can overcome the magnetic 
field of a star and the disc will revert to a Shakura-Sunyaev 
form. 

The remaining parts of the solution are also easily found 
and are collected below. In all cases they differ from the 
usual Shakura-Sunyaev form by the presence of the factor 



H = 1.7 x 10V 



-1/10^3/20 



M- 3/S R 9/S f 3 / 5 ,, 3 / 20 
16 IV1 1 n W J u 



p= 3.7xI0~ B a 
T c = 2.5 x I0 4 a 



-7/10 



-l/ 5M 3 /10 M y4^-3/4 /6 /5 ? 3/10 



ion -4/5 ,>l/5 ,4/5 1/5 
1 14 ^, 4 /5 A >3/10 



v = 2.9x 10 i V /o M 1 J 6 /lu M 1 - 1/ * J R u V7 



1/4 p 3/4 ,6/5 3/10, 



(37) 
u cm" 3 (38) 
K (39) 
(40) 
(41) 



The radial velocity of the disc material can be calculated 
directly from 



f R ~ 


3 
2 




Rt~ 


(2-7)\ 




.Rco - 






R. 







-%Lr-^ 

2tt 



(42) 



/ = 



(7-2) 
(7-1) 

R 



1 - 



Rt~ 




.R. 





1/2 



1/4 



Till 



(31) 



(32) 



(33) 



MR'' (7- 2) / 4 ' 
where i?* replaces R t throughout if R t < R±. In the limit 
= 0, this therefore yields the non-magnetic result in equa- 



2.5 The general opacity case 

Kramer's opacity is valid for discs where T c > 1 x 10 4 K. 
For discs around young stars, for example, this condition 
will not always be valid. The range of opacities in such discs 
is discussed in ISemenov et, al.l l)2003h . It is therefore useful 
to examine some alternative opacity laws and apply them 
in a manner similar to that found above. The process is 
straightforward, providing the opacity can be approximated 
by the form 
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t = Ko Ep a T c b , (43) 

where kq, a and b are constants. A collection of opacity 
prescriptions in t his form can be found in the appendix of 
iBell fc Lid lll994l) . The equations of the full solution are col- 
lected below, where d = 20 + 6a — 46. 




3 FORM OF FUNCTIONS 

Some typical disc structures will now be illustrated, using 
the full disc solution quoted in equations Ij44^ - Ij49^ . The ef- 
fect upon the solution of altering some of the parameters will 
also be shown. In every case the truncation radius, which is 
required for the solution, is obtained from equation (J2HJ- It is 
therefore implicitly assumed that 7 = 7/2. Fig.|2|shows the 
solution as applied to a disc surrounding a typical young 
stellar object (YSO). The disc is plotted from R = to 
R — 3.5 x 10 12 cm. This is an arbitrary maximum radius, 
used to illustrate the behaviour of the inner disc since, in 
a real YSO, the disc would be much larger than this. In 
the case of a single YSO, which occupies a simple potential, 
the model can be extended outwards indefinitely. It should 
be noted however that, in a model with a great radial ex- 
tent, it is unlikely that the entire disc co uld be correctl y 
represented with a single a parameter (e.g. lGammielll996l) . 



and a single opacity prescription. A further complication in 
these massive, relatively cool discs, is that the gravitational 
instability may also contribute towards the disc viscosity. 
In the present case however, where only the inner disc is 
modelled, Kramer's opacity is assumed throughout the sim- 
ulated region and a global value of a — 0.01 is adopted. 
This is undoubtedly a simplification, but is adequate for the 
purpose of illustration. This value of a in YSOs is consis- 
tent with FU Orionis r ecurrence times in the magnetic case. 
jMatthews et alJl20o3) . 

The sample star has typical properties for a low mass 
YSO so that i?* = 1 x 10 11 cm, K'h = 1 M , P 8pin = 3 d. The 
vertical magnetic field at the surface is B z — 500 G, which 
is quite low when compared to recent measurement s which 
give fields of order kG (e.g. I.Tohns-Krull et, al.Nl9 99'l. How- 
ever there is an uncertainty in how /3 relates to field strength, 
which makes comparison with measured field strengths im- 
precise. The disc is illustrated at mass transfer rates of be- 
tween M = 5 x 10~ 9 M yr" 1 and M = 2 x 10~ 8 M w yr" 1 
This is a typical M as measured at the star dGullbring et alj 
1998) and must, in steady state, be the M into every anmihis 
of the disc. 

Fig. [5^ shows surface density as a function of radius. 
The form is similar to that of the usual non-magnetic 
Shakura-Sunyaev solution. It differs however in that the in- 
ner disc is truncated further out than the stellar surface. In 
this plot only results at radii greater than R t are shown since 
within R t the solution becomes unphysical. In reality this re- 
gion may contain more complex accretion streams onto the 
magnetic poles of the star, which are difficult to reproduce 
in a one-dimensional treatment. Truncation should occur in 
the region where the magnetic time-scale and viscous time- 
scale become comparable. In this region the magnetic field 
is also sufficien tly influential to produ ce effects such as mag- 
netic warping dO'Sullivan et al1l2004l) . As the accretion rate 
increases, so does the surface density. As expected, an in- 
creased accretion rate also pushes the truncation radii in to 
smaller radii. 

The disc scaleheight (Fig. |^p) follows an approximately 
linear relation with radius, as with the usual Shakura- 
Sunyaev solution. However, truncation occurs some distance 
from the star in the same manner and at the same radii as 
with the surface density. The disc becomes thicker at higher 
accretion rates, but for realistic parameters remains of the 
order H ~ 0.1J? so that the thin disc approximation appears 
to be justified with these parameters. The temperature of 
the disc mid-plane, as illustrated by Fig.|5Ji, follows a series 
of curves similar in form to those of the surface density. This 
is to be expected since T c tx E 2 for Kramer's opacity. Fig. 

plots |i>r/ivI as a function of radius where v v is assumed 
to be Keplerian and «r is calculated using equation l|42|l . It 
can be seen that ur <C v v outside the truncation radius, so 
the assumption that the disc is cool is valid for the region 
which is modelled. 

The advantage of the solution derived in Section f2.5l is 
that it can be applied to a general opacity case, providing 
that the opacity law follows the relation given by equation 
l|43> . The discs about young stars are likely to be dusty 
and cooler than those in most cataclysmic variables or X- 
ray binaries. The surface density and temperature of such 
protoplanetary discs are plotted in Fig. [2] for three different 
opacity laws. The opacity prescriptions for those plots were 
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Figure 2. Analytic model of a typical YSO disc showing the surface density, disc scalehcight, viscosity and central temperature as a 
function of radius. The ratio of radial to azimuthal velocity is also shown. The disc is plotted from R = to R = 3.5 X 10 12 cm. The 
star has a radius R+ = 1 X 10 11 cm, mass M* = 1 Mq, spin period P sp i n = 3 d and a surface magnetic field of B = 500 G. The mass 
transfer rate throughout the disc are from M = 5 X 10 — 9 Mq yr -1 to 2 X 10 — 8 Mq yr" 1 . An alpha parameter of a = 0.01 is adopted, 
and Kramer's opacity is assumed. The vertical lines illustrate the stellar radius and the corotation radius. 
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Figure 3. Analytic model of a typical YSO disc showing surface density and temperature for three different opacity prescriptions. The 
disc is plotted from R = to R = 3.5 X 10 12 cm. The star has a radius i?* = 1 X 10 11 cm, mass M+ = 1 Mq, spin period P sp j n = 3 d and 
a surface magnetic field of B = 500 G. The disc has a mass transfer rate of M = 1 X 10 — 8 M^ yr -1 , an alpha par ameter of a = 0.01. In 
addition to Kramer's opacity (kq = 5 x 10 24 , a = 1, b = —3.5), two opacity prescriptions from lBell Hr, Linl ll994l) are illustrated: metal 
grains (fto = 0.1, a = 0, b = 0.5), and ice grains (fto = 2 X 10 — 4 , a = 0, 6 = 2). The vertical lines illustrate the stellar radius and the 
corotation radius. 




Figure 4. Analytic model of a typical dwarf nova disc showing surface density and viscosity for two different values of a and of M . The 
disc is plotted from i? = 0toi? = 2x 10 10 cm. The star has a radius = 1 X 10 9 cm, mass M* = 1 Mq , spin period P sp m = 100 s 
and a surface magnetic field of B = 2 X 10 3 G. Kramer's opacity is assumed. The alpha parameter takes values of a = 0.01 and 0.1, 
representing the disc in quiescence and outburst respectively. The disc has a mass transfer rate of M = 1 X 10 — 11 Mq yr _1 in quiescence 
and 1 X 10 — 10 Mq yr — 1 in outburst. The vertical lines illustrate the stellar radius and the corotation radius. 



obtained from lBell fc LinHl994ll . The shape of the curves are 
affected only slightly by the change in opacity, but the mag- 
nitude of the densities and of the temperatures are altered 
by several orders. It is notable, although unsurprising given 
the form of equation 12(il . that the truncation radius does 
not vary with opacity. Plotting the temperature of the disc 
mid-plane provides a useful check of consistency for opacity 
assumptions. It is clear for example that the inner disc here 
should not be represented by an opacity law which models 
the behaviour of ice crystals, since T c > 1 x 10 4 K in much 
of the inner disc if such an assumption is made. 

In addition to the above-mentioned uncertainties re- 



garding the disc viscosity, the young stellar objects have lim- 
itations as test cases for the new disc solution. The study of 
outburst behaviour is particularly difficult with YSOs. Since 
only the very inner part of the disc is modelled, the mass 
transfer rate into the simulated region would in reality be 
related to the viscous state of the outer disc so that a and M 
are not independent. This complication is avoided for dwarf 
novae, where the whole disc can be simulated. In this case 
the mass transfer rate into the accretion disc is unrelated 
to the viscous state of the disc. Moreover, the behaviour of 
dwarf nova discs is in general better understood than that of 
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Figure 5. Analytic model of a typical dwarf nova disc showing surface density and v as functions of radius for a range of magnetic field 
strengths. The disc is plotted from i? = 0toij = 2x 10 10 cm. The star has a radius i?* = 1 X 10 9 cm, mass M* = 1 Mq and spin period 
^spin = 100 s. The disc has a mass transfer rate of M = 1 X 10 Mq yr and Kramer's opacity is assumed. The alpha parameter takes 
the value of a = 0.01. The stellar magnetic field ranges from B = 2 X 10 4 to 1.15 X 10 4 G in the upper two plots and from B = 2 X 10 3 
to G in the lower plots. These values are selected so that the ratio /3/M is the same as in the three graphs plotted in Fig. [5] in order 
that the effects of M and /3 may be seen independently. 



those around YSOs l)Warnerlll995h . In the hot viscous state 
dwarf nova discs are very close to the steady state. 

In the case of binaries however, the disc occupies a more 
complicated potential and an axisymetric model such as the 
one described in Section [5] must, unavoidably, be an incom- 
plete one. Effects such as tidal resonances cannot easily be 
modelled in one dimension. A one-dimensional model should 
however be a reasonable approximation whilst magnetic and 
viscous forces dominate over tidal ones as they do in the in- 
ner disc. 

Fig.|I]shows the steady-state density profiles and viscos- 
ity profiles of a typical dwarf nova disc in a hot, high alpha; 
and in a cold, low alpha viscous state. The disc is plotted 
from R = to 2 x 10 10 cm while the central white dwarf has 
a radius of R* = 1 X 10 9 cm. The white dwarf is spinning 
with a period of P sp i n = 100 s and the surface magnetic field 
is B — 2 x 10 3 G, which would normally be insufficient for 
the binary to be considered magnet ic, since interm ediate po- 
lar binaries have fields of B ~ MG JWarnerlll995t) . It should 



be noted however that magnetic field strengths for white 
dwarfs with weaker fields are not generally well known. The 
disc opacity is modelled using Kramer's law throughout this 
and subsequent disc models. 

It is interesting, and possibly counterintuitive, that a 
change in a does not, in itself, cause R t to migrate. That this 
must be the case is very clear however from equation 12611 . In 
fact, for a given steady-state solution, the only parameters 
which set the position of the truncation radius are M, j3, 
R* and P sp m- However, the truncation radius is not truly 
independent of a since, during the change of the viscous 
state, M also changes and R t is a function of M. 

Although in reality the disc of a dwarf nova does not 
occupy a true steady state during either outburst or qui- 
escence, it is a reasonable approximation to represent these 
two phases by two steady states. The outburst cycle can then 
be interpreted as a cycle between these two steady states. 
The first plot in Fig. 0] illustrates the density profile of the 
disc in the hot and cold states where a = 0.1 and 0.01 re- 
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Figure 6. Analytic models of a typical dwarf nova disc with a range of mass transfer rates. The disc is plotted from R = to 
R = 2 X 10 10 cm. The star has a radius -R* = 1 X 10 9 cm, mass = 1 Mq, spin period P sp i n = 100 s and a surface magnetic field 
of B = 2 X 10 4 G in the first plot and B = 2 X 10 3 G in the second. The disc has a range of mass transfer rates from M = 1 X 10 
to 3 X 10 -11 Mq yr _1 and Kramer's opacity is assumed. The alpha parameter takes the value of a = 0.01. Again, the vertical lines 
illustrate the stellar radius and the corotation radius. 
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Figure 7. Analytic model of an X-ray binary disc containing a magnetic neutron star which rotates at different speeds. The disc is 
plotted from R = to R, = 8 X 10 8 cm, showing surface density and h as functions of radius. The star has a radius i?* = 1 X 10 6 cm, 
mass Mi, = 1.4 Mq and spin periods ranging from P S pin = 0.1s to 10s. The disc has a mass transfer rate of M = 1 X 10 — 9 Mq yr — 1 and 
Kramer's opacity is assumed. The alpha parameter takes the value of a = 0.01. The stellar magnetic field is set to B = 2 X 10 s G. The 
vertical lines, which are plotted with the same line style as their associated plots, represent the corotation radii of the primary, at the 
three different spin periods. 



spectively while M = 1 X 10 11 and 1 x 10 10 . The second 
plot shows disc viscosity for the same two steady states. If 
the viscous time-scale is given by 

ivisc = — (50) 

then in the hot state the viscous time-scale is reduced from 
ivisc ~ 1000 d to ~ 100 d at R = 2 x 10 10 cm. This allows 
viscous processes to remove the inner disc material more 
rapidly and reduce the build up of mass. 

It was noted in Section 12.41 that the new magnetic so- 
lution modifies the form of the standard Shakura-Sunyaev 
model only by the function v which in turn varies from unity 



according to the ratio /3/M, and by replacing R+ by R t . 
This means that when /3 tends to zero, and, Rt tends to 
i?*, we recover the Shakura-Sunyaev solution and that for 
high M the magnetic effect becomes less important. How- 
ever, since M also appears elsewhere in the solution, the 
effects of varying M and /3 are best examined separately. 
Figs. |3 and 1^1 show a cold dwarf nova disc, where the white 
dwarf has the same properties as in Fig. 2] In Fig. [S] the 
effect on surface density of a change in B is illustrated. In 
and^ B is varied from 2 x 10 3 G to 1.15 x 10 3 G. 
The function v is plotted for the same parameters in Fig. 
03d. As has been mentioned, there is no steady state for a 
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magnetic propeller. However the regime illustrated here is 
that of a 'near-propeller'. The truncation radius is very close 
to its maximum radius, that of the corotation radius, and 
moves very little with increased field strength. Therefore the 
great majority of the disc mass lies outside corotation and is 
propelled outwards. The stronger the field, the greater this 
tendency will be. If the magnetic advection term is to be 
balanced by a viscous diffusion term then, for an increasing 
magnetic field, an increased surface density is required. In 
this case the function v acts to enhance the surface density 
over most of the disc. Mathematically, there is no limit to 
the surface density enhancements which may be obtained in 
this manner. Physically however, there are two reasons why 
this is not so. Either the disc will, at some point, reach so 
high a surface density that a thermal- viscous outburst will 
occur, or the disc will become very thick and the thin disc 
treatment will cease to be valid. 

For lower field strengths the disc behaves differently. 
This 'strong accretor' regime is illustrated in Figs. and 
QJl. Where (3 = we have the Shakura-Sunyaev solution and 
v retains a value of unity at all radii. As j3 is increased the 
truncation radius migrates outwards and in this case the 
surface density decreases throughout the disc. The effect is 
however more pronounced in the inner disc and the solution 
tends to the non-magnetic case with increasing R, as can 
clearly be seen from the form of v, which never exceeds unity, 
and so always reduces E. 

Fig. H]] shows how the steady-state solution for a dwarf 
nova system, which is otherwise identical to that shown in 
Fig. varies when M is changed and (3 is held constant. 
In both plots M is varied from I x 10~ n M Q yr _1 to 3 x 
10 -11 Mq yr _1 . Care is taken that the ratios f3/M are, in 
the first plot, the same as in Fig.|^, and, in the second plot, 
are the same as in Fig. 0:, although the 5 = plot has no 
analogue here. It can be seen that, while the truncation radii 
are identical to those in the case where B is decreased, the 
surface densities are higher in the case where M is increased. 
In fact even in the 'strong accretor' case, surface density 
increases with increasing M . This is as expected from the 
form of equation (|44pi . It is significant that for M > 2 x 
ICP 11 Mq yr _1 the truncation radii obtained by magnetic 
truncation a re comparab le to those of order R t ~ I0 9 cm 
predicted bv lKina 1(1997) as a result of irradiation. However 
at lower mass transfer rates the magnetic R t will become 
larger than that due to irradiation. 

In order to illustrate the effect of a change in P sp in 
upon the disc, a more rapidly rotating star is advantageous. 
For this reason we use a neutron star in a low-mass X-ray 
binary disc. In Fig. the disc is plotted from R = to 
R = 8 x I0 8 cm. This is only a model of the centre of the 
disc, since this is where the solution diverges most from the 
non-magnetic case. The central star has typical properties 
for a neutron star in a low mass X-ray binary, with radius 
7?* = I x 10 cm, mass = 1.4 Mq and a range of spin 
periods from P sp m = 10 s to 0.1 s. The stellar magnetic field 
is set to B — 2 x 10 8 G for all plots. The mass transfer rate 
has a constant value of M = 1 x 10~ 9 Mq yr _1 . The disc is 
assumed to be in a cool state and to obey Kramer's opacity 
throughout. Since we only model the centre of the disc it is 
reasonable to ignore tidal effects. The corotation radii here 
vary according to P sp i n and are marked as vertical lines with 
the same line style as their associated plots. 



The surface density plot shows that an increased spin 
rate causes the magnitude to increase and brings R t closer to 
the star. To understand what causes this it is useful to exam- 
ine the form of h, defined in equation (I31L which represents 
the potential for departure from the non-magnetic model as 
a function of R and of R co . This function is independent of 
the (3/M ratio, which is in any case held constant through- 
out Fig.0 Thefunct ion is made up of two terms, the first of 
which depends upon the stellar spin and the second of which 
does not. The two terms have opposing signs for 7 > 2. The 
expression can be interpreted so that the second term repre- 
sents the residual accretion effect of a notional non-rotating 
magnetic star. For long spin periods, and hence for large 
values of P co , it is clear that h tends to this non-rotational 
case. The first term in h is the propeller-like term which 
becomes more important when R > R co . The transition be- 
tween accretor and propeller is not however so abrupt or 
straightforward as might be imagined because of the differ- 
ent powers of R in the two terms. The first case in Fig. |7| 
represents a 'near-propeller' solution, which has the signa- 
ture of a negative h, which increases surface density. In this 
case a 'near-propeller' solution requires a stronger field than 
a 'strong accretor' solution to truncate the disc at the same 
radius. This may seem counterintu itive, but it agrees with 
the results of lLivio fc Pringlel il992T) and can be understood 
by comparing the mechanisms of the magnetic propeller and 
magnetic accretor. With a given spin period, and disc mass, 
the magnetic propeller can potentially clear a larger hole in 
the disc, because it is not limited by the corotation radius. 
However, because a magnetic propeller has to act in oppo- 
sition to viscous processes, and because it causes a build up 
of mass at the edge of the hole, a stronger field is required 
to maintain the same size of hole in the propeller case. For 
faster spinning stars, the corotation radius is smaller and so 
the propeller mechanism operates from closer to the star. 
This causes the truncation radius to become smaller. 



4 DISCUSSION 

A new self consistent analytic solution has been developed 
for a thin accretion disc, under the influence of a central 
magnetic field. Since this model represents a steady-state 
disc it is not applicable to all accretion discs, and cannot 
be used for the magnetic propeller. The model was obtained 
in a similar manner to the non-magnetic Shakura-Sunyaev 
model, and indeed the result tends to that non-magnetic 
case when the magnetic parameter (3 tends to zero. The 
two models also converge at large radii. The analytic model 
has also been expanded to include the effects of a range of 
opacities in addition to Kramer's law. The solution has been 
applied in the cases of young stellar objects, cataclysmic 
variables and X-ray binaries. 

The new model produces a solution in which the in- 
ner disc is frequently truncated further out than the stellar 
boundary layer. This is caused by magnetically induced ac- 
cretion. The results depend upon how the position of the 
truncation radius -Rt varies as a function of the parame- 
ters of disc and star. It has been shown that the trunca- 
tion radius is independent of the opacity prescription used. 
The magnetic interaction model is uncertain, so that nei- 
ther the exact values of the magnetic parameters f3 and 7 
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nor their dependence on the field strength of the star are so 
far well known. Further work is therefore required to better 
constrain the relationship between the surface field and the 
truncation radius. In addition, a comparison of the position 
of the truncation radius between the one dimensional so- 
lution and three dimensional techniques, such as smoothed 
particle hydrodynamics, could be used to ensure that the 
one dimensional inner boundary condition acts in the same 
manner as the three dimensional analogue. 

The critical radius at which thermal viscous 
outbursts begin varies almost linearly with radius 
flCannizzo. Shatter fc Wheelerl Il98ftf) . Therefore a trun- 
cated disc can store more mass before it goes into outburst. 
As a result magnetic fields can exert a strong influence 
on outburst cycles in accretion discs by forcing outbursts 
to begin furthe r out in the d is c. Th is is discussed in 
more detail by iMatthews et alJ J2004J) . It would be in- 
structive to perform full outburst simulations in which 
disc annuli are permitted to switch between the hot and 
cold state according to local triggers, and on a thermal 
time-scale. Similar work has been carried out to model 
dwarf nova outbursts in smoothed par ticle hydrodynamics 
(e.g. iTruss. Wvnn fc Wheatlevl l2004h . The steady-state 
viscosity prescription should not be applied to a time 
evolution model because it becomes inaccurate at low 
densities. In addition a negative value of v, which can occur 
for R < R t , will cause equation l)14|l to have sinusoidal 
solutions. More realistic simulations could be performed 
by relating viscosity to surface density or by using a more 
physically motivated viscosity such as t hat generated by 
the m agneto-rotational instability (e.g. iBalbus fc Hawlevl 
Il99ll) . 

Additional future work will involve the generalisation of 
the new solution to encompass torques from other sources 
such as planets embedded in the disc. This would have im- 
portant implications for the study of planet formation and 
migration in discs. The effect of the magnetic torque on the 
central star's spin could also be investigated in detail. The 
long term effect of such spin evolution may be found to have 
an significant effect on the structure and behaviour of the 
disc. 
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